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We study scalar-tensor theory, k-essence and modified gravity witii Lagrange multiplier constraint 
which role is to reduce the number of degrees of freedom. Dark Energy cosmology of different 
types (ACDM, unified inflation with DE, smooth non-phantom/phantom transition epoch) is re- 
' constructed in such models. It is demonstrated that presence of Lagrange multiplier simplifies the 

' reconstruction scenario. It is shown that mathematical equivalence between scalar theory and F{R) 

, gravity is broken due to presence of constraint. The cosmological evolution is defined by the sec- 

■ ond F2{R) function dictated by the constraint. The convenient F{R) gravity sector is relevant for 

local tests. This opens the possibility to make originally non-realistic theory to be viable by adding 
the corresponding constraint. A general discussion on the role of Lagrange multipliers to make 
higher-derivative gravity canonical is developed. 
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I. INTRODUCTION 

The understanding of late-time cosmic acceleration (the so-called Dark Energy (DE) epoch) is one of the main 
^ ■ challenges of modern cosmology. There is no yet totally convincing theoretical DE model which is related also with 
, the lack of precise observational data. This requests the introduction of new DEs in order to present fundamental and 
simple DE. It seems that standard strategy of combination of known models does not work and some new approach 
should be invented. Recently, a quite interesting new DE model P, [3| was proposed which consists of two scalar fields 
^ I where one of scalars represents the Lagrange multiplier. This multiplier puts the constraint of special form on the 
^— I ' second scalar field. As a result, the whole system contains the single dynamical degree of freedom. It was shown that 
0^ 1 energy of such theory flows along time-likegeodesics which is similar to the dust. Nevertheless, the theory contains 
■ non-zero pressure (dusty dark energy) [l], [J. It may be very natural candidate for unification of Dark Energy and 
\ Dark Matter. It is interesting to study the role of such novel construction in the known DE models because it may 
' completely change the structure of theory, its cosmological solutions and cosmological perturbations theory and may 
give better and/or qualitatively different fit towards to observational data. Moreover, the role of Lagrange multipliers 
may be twofold because they are often used for reducing of higher- derivative gravity systems in canonical formulation. 
" , [ ' Hence, the most natural area for the study of Lagrange multiplier constraint is modified gravity. In addition, adding 
^ , such constraint in modified gravity may significally improve the ultraviolet properties of the covariant theory [3| 

leading to its renormalizability conjecture. 
^ In the present letter we discuss the role of such Lagrange multiplier constraint in modified gravity. This is the 

first study of this sort, so the number of models are discussed in order to understand the general role of Lagrange 
multiplier construction and its impact to cosmology. We start from unified (phantom/canonical) scalar theory with 
Lagrange multiplier constraint. The reconstruction scenario for such theory as well as for k-essence theory is developed 
and shown to be significally simpler than in the case without constraint. Dark Energy cosmology of different types 
is reconstructed including standard ACDM, unified inflation-Dark Energy era and non-phantom/phantom transition 
cosmology which turns out to be smooth. The comparison with standard theory is done, it is shown that different 
potentials describe such cosmological evolution. Moreover, convenient scalar does not propagate and does not generate 
the extra force. Despite the fact that theory is described by single degree of freedom, it may be represented in BD-form 
where Lagrange multiplier Iras ghost-like kinetic term. 

The relation of such theory with modified F(i?)-gravity is investigated. It is shown that mathematical equivalence 
between scalar theory and F(i?)-gravity is broken because of the presence of Lagrange multiplier constraint. New 
constrained form of modified gravity is introduced. It consists of standard i^i (i?)-term and Lagrange multiplier which 
constrains kinetic- like term for curvature by second function F2{R). The cosmological dynamics of such modified 
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gravity depends only from the choice of F2{R) function. Dark energy cosmology may be easily reconstructed using 
its form. 

Unlike to scalar theory, the Lagrange multiplier propagates. Hence, the role of Fi{R) is in modification of Newton 
law which may cause the extra modification of constraint in order to reproduce the standard Newton law. This is also 
novel property, as in convenient modified gravity if theory does not pass newton law regime it is ruled out as realistic 
one. 

Another application of Lagrange multiplier in canonical formulation is discussed. General method of reduction of 
higher derivatives is developed. It can be adopted for any higher-order theories. Specifically, gravity theories of order 
{2k + 4), with k an integer, need (fc + 1) Lagrange multipliers to be made canonical, as we will discuss below. 

The layout of this letter is the following. In Sec. II, we study scalar-tensor and k-essence theories with Lagrange 
multiplier constraint. DE cosmologies which are generated by the self- interaction potential and Lagrange multiplier 
constrained fiuid are investigated. Sec. Ill is devoted to the same issue but F{R) gravity is considered. In this case, 
cosmological dynamics is qualitatively changed due to lack of mathematical equivalence with scalar-tensor theory and 
presence of second F2{R) function caused by Lagrange multiplier constraint. In Sec. IV, general considerations on 
Lagrange multipliers are developed. We show that a given gravity theory of order (2fc + 4) can be made canonical by 
introducing [k -\- 1) suitable Lagrange multipliers. Conclusions are drawn in Sec.V. 



II. DARK ENERGY IN SCALAR THEORY WITH LAGRANGE MULTIPLIER 

In this section, we consider the accelerating FRW cosmology in the theory with two scalars where one of scalars is 
Lagrange multiplier which constrains the field equation of second scalar. The starting action has the following form: 

j d'x^g I A _ ^d^<t>d'^^ - V{4>) - A l^-d^4>d^^<t> + U{4>)^ I . (1) 

Here A is the Lagrange multiplier field. Depending on the sign of the potential w(0) (see ref.[3|) the first scalar could 
be the canonical scalar or the phantom. The Einstein equation has the following form: 

^ {R,.^-\g,uR) - \g,u j-^gp^y^ - v{<j>) - A (\dp<j^d^<j> + u{^)] ] + ''^'^1,+ ^ d.^d^. (2) 



We now work in the FRW metric with fiat spatial part: 

ds^ = -dt^ + a{tf {dx'f . (3) 



i=l,2,3 

Then by the variation over A, we obtain 

i2 



= y - U{<P) . (4) 



The FRW equations are given by 



3_^2 ^ + ^ y^^^ ^ ^f^(^) ^ (^(^) ^ 2A) U{cf>) + Vicj,) , 

These equations show that the EoS-parameter has the following form: 



Eq.Q can be integrated as 



""^^ [u{ct)) + 2A) U{(t)) + V{(t)) ■ ^ ' 



as long as [/(</>) is positive. Then by solving Eq.Q with respect to 4>, one can find the t-dependence of ^: (j> = 4>Q{t). 
Substituting the expression of (l)o{t) into ([5]), we obtain a differential equation for H, which gives the ^-dependence of 
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H: H = Ho{t). Finally substituting (f'oit) and Holt) into we can find the t-dependence of the Lagrange multiplier 
field A: 

^i<f>o{t)) , 1 f 3 „ , ,2 



Note that due to the constraint (|4]), the dynamics is completely changed. Usually if the scalar field exists, the 
propagation of the field generates an extra force and often violates the observational constraints on the Newton law. 
Due to the constraint (jlj, however, there does not appear the propagating mode of 0. To see this, we consider the 
perturbation from the solution (j)o{t) given by (O, 

4>^Mt) + H- (9) 

Here we assume that 64> can depend on both of the time coordinate t and the spatial coordinates. Then the constraint 
equation ([4]) gives, 

dt ^ ^ ^ ^ 

dt 

When U'{(j))/^^ > 0, the perturbation Scjj grows with time and therefore the solution (/jq becomes unstable. On the 
other hand, when U'{(j))/^^ < 0, the solution 0o becomes stable. Anyway there does not appear the oscillating mode 
in (50 and therefore Scj) does not propagate nor does not generate the extra force. We should also note that in the 
equations ^ and ([SJ, there does not appear the term containing the derivatives of A, which shows that the multiplier 
field does not propagate. We now consider the reconstruction, that is, when the behavior of the Hubble rate H: 
H = Ho{t) is known, how can we construct the action in the form of ([T]) to reproduce Ho{t). First we choose U{(j)) 
appropriately so that we can easily integrate Eq. ([7]) to find explicit form of 0o(t) and t = t{(j)). Then Eq. ([5]) gives 

T/(0) = ^ (2Ho (t (0)) + 3Ho [t (</>))') + ^(0)C/(0) . (11) 

Here uj{ip) can be arbitrary. Hence, one finds that choosing V{(l)) as ([TT|) . the arbitrary Hubble rate H = iJo(i) can 
be reproduced. 

Some examples may illustrate the reconstruction scheme. Just for simplicity, we assume 

4 

^(</>) = l, f/(0) = — . (12) 

Here to is a constant with the dimension of mass and canonical scalar is considered. Then Eq. ([7]) tells 

cl) = m^t. (13) 
Here we have chosen the + sign in ± in ([7]). Then Eq. PT|) has the following form: 

Eq. ([5]) indicates that 

X = -l-^Ho(^)- (15) 

\ TO^ / 

Let us consider Ho{t) corresponding to the ACDM model: 

i/o(0 = |cothQ . (16) 

Here I is the length parameter related with the cosmological constant. The potential V{4>) becomes a constant 

4 
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which may be regarded as a cosmological constant. The scalar field plays the role of dust[l|, 0|. Thus, one obtains 
dark energy produced by dusty Lagrange multiplier field. Following the proposal of ref.pj, one can reconstruct the 
model unifying the inflation and the late-time acceleration: 



72- 



Ho{t) = ^ ° . (18) 

When t 0, H behaves as H = Hj + O (t^)- Then the universe is almost de Sitter space corresponding to the 
inflation. On the other hand, when t — > oo, H behaves a.s H = Hl + O and the universe becomes asymptotically 

de Sitter space describing dark energy era. Applying the above procedure one gets: 

^ ^ 'f^^. ^ + ^ . (19) 

As a final example, we consider the dark energy model which admits the transition [5] from non-phantom phase to 
phantom phase: 

Hoit) ^ho(l + . (20) 



Then we find 



y,^. ^ I {3hl - 2ho) {3hl - 2ho) Gfe^m^ \ 

The theory with above scalar potential describes the transition from non-phantom phase to phantom era. The EoS 
parameter 

^i'0A--l-^rT— ^, (22) 

which shows the phantom crossing (w — ~1 crossing) when t — ts/2, that is, w > —1 when t < ts/2 and w < —I when 
t > ts/2. Eq. (|20p demonstrates that there is a Big Rip singularity at t = tg (for first works on Big Rip singularity, 
see 0). To understand better the role of Lagrange multiplier constraint we consider single scalar-tensor theory 



^d.cpd'^cj^-VicP)} . (23) 



The reconstruction can be performed by choosing 



w(0) = -^i?^(0) , V{cb) = 4 i^Hoi^f + H'M)) ■ (24) 



Then it follows 



d^ = t, H^f{t). (25) 
Especially in case of non-phantom/phantom transition (PO)) . one gets[l| 

The form of the potential is similar to that in (PT|) but the coefficients are different from each other. We should also 
note that in (PT|) . the mass dimension of the scalar field is unity as standard but that of the scalar field in is 
minus unity by following fl, 0|- The resemblance of two potentials pT|l and ((26|) comes from the expressions ([Tl]) or 
more explicitly (|14p and ([M)) . Note that in (HJ]), the scalar field is chosen, essentially to be time coordinate as in 
by the choice of U{4>) in (fT^ . Expression pTjl comes from the effective pressure but comes from the difference 
between the effective pressure and the effective energy density. Then the coefficients in two potentials are different 
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from each other. In spite of the similarity in the two models (|23p and there is a big difference in them, that is, 
the scalar in ([2^ propagates and often gives a correction to the Newton law but the scalars in ([T]) do not propagate 
and do not give any correction to the Newton law. Now both of the scalar-tensor theories with a constraint ((2T|) 
and without a constraint (j26p give the identical EoS parameter w: w — Wipx in (|22p . Then in both of the models, 
the phantom crossing occurs when t = is/2 and there is a Big Rip singularity dX t — ts- We should note that uj{(f)) 
vanishes when the phantom crossing occurs at i = ts/2, which shows that if we redefine the scalar field as 



ip = j d4>^\ujm , (27) 

the action can be rewritten as 

j d'x^g I T ^S^-^^V ~ V{^) I . (28) 

The sign in front of the kinetic term depends on the sign of uj{(j)). Therefore at the point of the phantom transition, 
the sign changes. In this sense, in the model (|23p . the phantom transition is not smooth and it has been shown that 
the transition is very unstable ^ Jj . On the other hand, for the model with a constraint ([I} , the transition seems to 
be smooth. Finally in this section, we show that the action ([IJ can be rewritten as the Brans-Dicke theory coupled 
with a (phantom) scalar field. Let Lo{(j) — 1 and consider the conformal transformation 

5m- = e'^ , e'^'^ = (1 + A)^' . (29) 
Then the action ([1]) has the following form: 

S = j d^x^ I {r + \d^^ad^a^ - \d^<pd^<j, - e^^/^y (<^) - (e^^/^ - e^/^) U{ct>) | . (30) 

Now the constraint disappears and the scalar field (j) ^'^d the ex-Lagrange multiplier field A or cr seems to be prop- 
agating. However, the kinetic term of a is not canonical (ghost-like). Then effectively the amplitude coming from 
the propagation of (p might be canceled by the propagation of a. One may choose V{(j)) so that the smooth phantom 
crossing occurs in our model. Thus, we demonstrated that the theory with two scalars where one of them is Lagrange 
multiplier may lead to variety of dark energy cosmologies. Despite the fact that the theory possesses single dynamical 
degree of freedom, its cosmology seems to be qualitatively different from the one of the theory with single scalar. 

It is well-known that scalar-tensor theory is mathematically equivalent to F{R) gravity which does not lead to 
physical equivalence of two theories^. Now we may try to transform the action ^ into the F(i?)-gravity form. 
Making the conformal transformation 

^ e^^Vf^^, , (31) 

the kinetic term of <b is canceled and one obtains 



S= f d^x^l '^ -e^"'^Viy(^)-A( ^e"^^9^(/.g^c/.-e^"^^[/(0)) I . (32) 




In the standard case, there is no last term in the above action and scalar becomes the auxiliary one (no the kinetic 
term). Since the term containing A includes the derivative of 0, we cannot integrate and/or delete the scalar field 4> 
and therefore it is difficult to rewrite the action ([T]) in the F(i?)-gravity form. Hence, the mathematical equivalence 
between two theories seems to be broken due to presence of the Lagrange multiplier term. 

As an extension of the scalar field with constraint, we may consider k-essence model[9] with a constraint: 



5 = y d4a;V^|^+X(0,X) + A(x-C/(0)) +L„,attcr|, X = <f>d ^cj^ , (33) 

where K is an appropriate function of 0. Let us check if the introduction of Lagrange multiplier in such theory does 
not destroy its consistency and the reconstruction discussed in [10]. The FRW equations are given by 

l.H^ = 2X^^^^-K{cj),X) + 2\X + p^,tt.r{t), --^(2H + 3H^^ ^ K {(f>, X) + p^,tter{t) . (34) 
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Here we include the contribution of matters with a constant EoS parameters Wi . Then the energy density Pmattcr and 
the pressure Pmattcr of matters are given by 

Pmatter = P^^^-^^^^^^^ , Pmattcr = ^^PO.a-'^'+"''^ ■ (35) 

i i 

Then the variation of ([55)) with respect to A gives, 

X - U{(j)) = . (36) 
We choose U{(j)) — m'^/2. So gives cj) = rn^t, then one can rewrite the equations in p4p in the fohowing form 

K {mh, mV2) = [2H + 3H^) - ^ w^po^a-^'^^+'"^'> , (37) 

dK (mH,X) 



= -^ij-A-^(l + w,)Po.a"'^'+"''^ (38) 

X=m4/2 i 



dX 

Then by using an appropriate function g{(f)/m'^), if we choose 

K{mh,my2) = (25"(^/m2) + 35'(</)/m2)') - ^ ^,po.ao"'^'+'"'^e-3(i+"'-)«(^/™^) , (39) 

i 

one gets the following solution of ([?7|) : 

H = g'{t) (a = aoe^(*)). (40) 

Note that X-dependence of K{(t), X) can be arbitrary as long as K{(t), X) satisfies ([5^ . The i-dependence of \{t) can 
be determined by using (|38p as 

-^i?-5](l + ^«0Po.a-3(i+-.). (41) 



dK{mH,X) 

\ — — 



dX 

In case of if (0, X) = K{X) action has the following form 



5 = y"d4xy^|^ + if(X) + A(x-C/((^)) + Wterj, X = -ia^09^</.. (42) 

If we set U{4>) = we find = 771^^2/2 and therefore U ^ X ^ Then the second equation of jMl) 

determines the form of K{X) as follows, 

K{X) = -\{2H(.j2Xj^) + 3i?2(y2XA^)) - u;,po^a(\/2XA^)-3(i+»') . (43) 

By differentiating the second equation of p4p . we find 

^tif'(m3iV2) ^-^{2H + 6HH) + 3i/^ii;,(l + w^)po^a-^^^+'^'^ . (44) 



Then combining (H^ with the first equation of (IMl) . the time dependence of A follows: 

A = - f ^ + V(l + w,)po,a-'^'+'^A -l-^+(^{H + mk) - 3il V «;,(1 + w,)po,a-^^^+'-A -\ . (45) 

Using above technique the specific examples of DE cosmology may be reconstructed. Theory remains to be consistent 
but the reconstruction examples qualitatively change. This means that the same potential which was used to produce 
given cosmology leads to different cosmology in the presence of Lagrange multiplier. Moreover, explicit realization of 
reconstruction scenario turns out to be significally simpler than without constraint. 
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III. FRW COSMOLOGY IN F(i?)-GRAVITY WITH LAGRANGE MULTIPLIER 

In this section we consider i^(i?)-gravity where Lagrange multipUer is introduced in the same way as in scalar theory 
of previous section. In the usual F(i?)-gravity, there appears the scalar mode called scalaron, which often affects the 
Newton law. In this section, we try to suppress the propagation of the scalaron by imposing the constraint under 
the Lagrange multiplier field. As a result, however, there seems to appear the propagating mode in the Lagrange 
multiplier field, which may break the Newton law but in somehow easier way. The solution of this question may 
request the additional modification of constraint. Another purpose of this section is the reconstruction. In usual 
F(i?)-gravity, we need to solve the complicated differential equation Q to realize the reconstruction program. In this 
section, we show that the reconstruction can be done more easily in the model with the Lagrange multiplier field. 
The starting action is given by 



S - / d^xV^ i.Fi{R)-X{ -d^Rdf'R + F2{R)]} . (46) 



Here A is the Lagrange multiplier field, again, which gives a constraint 

f 

On the other hand, by the variation of the metric g^^, we obtain an equation corresponding to the Einstein equation: 



-a^i?a^i? + F2(i?) = 0. (47) 



= ^gt^uFiiR) + ^d^Rd^R + {-R^, + V^V, - 5^,V2) {FUR) - XF.^{R) - (AV^i?)) . (48) 
If the Ricci curvature is covariantly constant and the scalar curvature is a constant: 

R^.. = ^df.. , R^Ro, (49) 

Eqs. (gT]) and (gS]) reduce to 

= F2(i?o), (50) 

= Fi{Ro) - ii?o (FiiRo) - XF^iRo)) . (51) 

If Eq. (f50|) has a solution, Eq. ([STI) can be solved with respect to the Lagrange multiplier field: 

, _ -Fi(i?o) + RqFI{Ro) 

Then if Rq is positive the above solution describes de Sitter space-time which may correspond to dark energy or 
inflationary epoch (for the proposal of gravitational unification of infiation with dark energy in modified gravity, see 
[ll|). For spatially-flat FRW metric Eqs. ((47|) and (^, ~ (0, 0)-component of p8)) have the following form: 

= -^R^ + F2{R), (53) 

= -\f,{R) + 18A [h + 4HHy + U (h + H^) 2,Hj\ {f[{R) XF^{R) + ("^ + 3ff] (x^) \ . (54) 



When F2{R) > 0, Eq. may be solved as 



which can be solved with respect to i? as a function oi t R = Fj^{t). Since 



I -^3^^ (55) 



dH ^ 

-R = 6— + 12i^^ (56) 



one can find the behavior of = ^ by solving the differential equation 

dH 



6^ + UH^ = Fnit) , (57) 
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By using the obtained solution for H = H{t) (and R = FR{t)), Eg. ([5^ becomes a differential equation for the 
multiplier field A and we can find the behavior of A, A = A(i). 

Conversely when the behavior of H(t) is known from the observational data, one may reconstruct F2 (R) to reproduce 
the behavior of H{t) by using ([55)) . H{t) gives the behavior of i? as i? = R{t), which can be solved with respect to t 
as t = t(R). Using ([55)l . the explicit form of F2{R) is found to be 



1 fdR"- ^ 



(58) 

t=t(fl.) 



Note that Fi (R) can be arbitrary function. Then the reconstruction of model can be more easily performed than that 
in the usual i^(i?)-gravity. As an explicit example, we may consider 

H{t) ^ ^ , (59) 



t 



where ho is a constant. Then 



-6/10 + m§ ^ -6ho + 12hl 
R= ^ or t=]J ^. (60) 



And therefore, we find 

dR 12 {-ho + 2hl) 2Ri 



dt ^6{-ho + 2hl) 

which gives 



(61) 



Another example is given by 



i? = — (1 - tanhwi) + ^ (1 + tanhwt) . (63) 



Here R± and uj are constants. Then t ±00, R — ±R± and therefore the space-time becomes asymptotically de 
Sitter. One may identify the epoch of i — !■ —00 as inflation and t +00 as late acceleration. Since 

2 cosh^ cot 2 y {R_-R+f J 

from Eq. (|58p . one gets 

{R^R+fu^f {R_+R+-2R)' , 

Thus, the unification of early-time inflation with dark energy epoch is possible also in constraint modified gravity. 
Hence, the universe evolution only depends on the constraint equation (|53p but does not depend on Fi{R). Fi{R) 
can only affect the correction to the Newton law. In convenient Fi{R) cosmology the whole dynamics is defined 
by the form of this function. With the constraint Fi{R) becomes irrelevant. The cosmological dynamics is 

defined by the form of F2{R). In _F(i?)-gravity, there appears the propagating mode, which is often called scalaron 
and which often violates the Newton law. In the same way as in the scalar-tensor theory around we may show 
that the scalaron does not propagate. In case of F(i?)-gravity, however, Eq. contains the second derivative of 
the multiplier field A although the Einstein equation ([2]) of the scalar theory with a constraint ([1]) does not contain 
the derivative of the multiplier field A and therefore A can be solved algebraically as in ([S]). In case of F(i?)-gravity, 
we need to solve the second order differential equation to find A, which might indicate that A could propagate and 
there might appear the correction to the Newton law. The magnitude of the correction could depend on the choice of 
Fi(i?) and/or i^2(-R)- 
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In order to investigate the Newton law, we choose Fi{R) as the Einstein one, 



F^{R)^:^, (66) 



and introduce the matter. Then Eq. (pS)) has the following form: 

= ^ Q.9M-i? - + ^T^. - + V, - V^) {XF^ (R) - (AV^i?)) . (67) 

For the solution where A = 0, Eg. pSj) reduces to the Einstein equation, 

1/^1 

k2 I 2' 



= 72 ( o5m-^ - ) + ■ (68) 



Here T^^ is the matter energy-momentum tensor. In the case without matter T^i, — 0, the Schwarzschild space-time, 
where R — R^^ = is a solution, which also satisfies the constraint equation (|47p if ^2(0) = 0. In case with matter 
T^v 7^ 0, however, the Einstein equation ([55)) gives 

R^-K^T. (69) 
Here T is the trace of the energy-momentum tensor. The constraint equation (|47p is rewritten to 

= y a^ra^T + f^ {-k^t) , (70) 



which is not always satisfied. Hence, in the presence of the matter, the constraint equation (|47p should be modified 
to be 

= ^df.Rd^'R + F2{R) - y^^Ta^T - F2 {-t^T) . (71) 
This indicates that the total constrained action with matter should be, instead of (^Sl 

S = j d^x^ ^-xi^d^Rd^'R + F2{R)-'^d^Td>'T-F2{-K'T)^+Cr^,,,,, , (72) 

In this case, the Newton law could be easily reproduced. Here >Cmattcr is the Lagrangian of the matter. Of course, 
qualitatively other form of constraint may also solve this problem. We also note that the form of the constraint could 
be correct when Fi{R) is given by For general Fi{R), the constraint could be changed. 
Since T vanishes in the vacuum as in the bulk of the universe, the constraint (|7ip reduces to 



= -df.Rd'^R + F2{R) ~ F2 (0) . (73) 



If F2 (0) — 0, the constraint ([7T|l gives (IT7|) and and the cosmological evolution could be generated. Note that 
F2 {R) in ([62]) satisfies the condition F2 (0) = but F2 {R) in ([65]) does not. In case the condition F2 (0) = is satisfied, 
there are two classes of solution in the constraint (|53p . One is a trivial solution i? = and another corresponds to 
the non-trivial cosmological evolution given by Near the solar systems and the galaxies, the solution could 

correspond to the trivial one i? = in order to reproduce the Newton law but in the bulk of the universe, the solution 
should correspond to ([55)) so that the evolution of the universe could be generated. It is not so trivial to show or to 
deny that the two solutions could be connected in the intermediate region between the region near the solar systems 
or galaxies and the region of the bulk universe. Maybe we need more careful (possibly numerical) analysis, which 
requests a future investigation in this direction. 

Thus, we demonstrated the possibility to describe the cosmological dynamics, including dark energy era, in modified 
-F(i?)-gravity with Lagrange constraint. It turns out that reconstruction which produces the viable cosmology in this 
case is realized qualitatively easier than in the convenient modified gravity. Moreover, to pass the local tests may 
request the additional modifications of constraint equations as is seen for emergence of Newton law regime. This may 
be caused by the fact that Lagrange multiplier propagates in such a theory. It is quite interesting observation because 
modifying the form of constraint one may arrive to qualitatively different predictions about local tests, which can 
make the same theory with different constraint to be viable! 
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IV. LAGRANGE MULTIPLIERS IN GENERIC HIGHER-ORDER THEORIES 



Let us discuss other application of Lagrange multiplies. In the optimization problems, such a method allows to 
find out extremal points (maxima and minima) of a function where one or more than one constraints are present. 
Specifically, Lagrange multipliers allow to calculate stationary points of the constrained function. In other words, the 
method reduces the search for stationary points of a n-variable function with /c-constraints to find out the stationary 
points of free (non-constrained) function of 71 + fc-variables: it introduces a new (unknown) scalar variable, the 
Lagrange multiplier, for each constraint present in the problem defining a new function (the Lagrangian) in terms 
of the original function, the constraints and the Lagrange multipliers. Up to this point, the Lagrange multipliers 
have been imposed a priori to modify the dynamics and select the form of the effective potential. Furthermore, by 
integrating the multipliers, cosmological solutions have been achieved. On the other hand, it is possible to show that 
the Lagrange multipliers are constraints capable of reducing the dynamics in higher order theories. Technically they 
are anholonomic constraints being time-dependent. They give rise to field equations which describe the dynamics of 
the further degrees of freedom coming from higher order theories. This fact is relevant to deal with such new degrees 
of freedom under the standard of effective scalar fields. 

With these considerations in mind, let us take into account generic higher-order theories described by the action 



A = 



^F{R, OR, a^R, tfi?) , 



(74) 



The field equations are 



R''" - l-g'^'R = i I ig''-(F - GR) + (5''^"" - g'^'^g^nO- 
z y [2 



where 



i=i j=i 

flU Act 

g g 



Act + 



, dF 
dUiR 



(□^■-l);,CT-^ 



dF 
dD'R 



(75) 



dF 

dam 



(76) 



These are pure gravity (2fc + 4)-order field equations. Matter can be taken into account by introducing, as above, the 
energy-momentum tensor T^i, . 

In order to better discuss the role of Lagrange multipliers, let us consider, for the moment, actions containing up 
to OR terms. In this case, we have eight-order field equations which becomes of sixth-order if the theories is linear 
in DR. If we take into account FRW point-like actions, we can reduce to the Lagrangian 



£ = £(a,a,i?, R, DR,{nR)), 



(77) 



by which one can deduce the Euler-Lagrange equations corresponding to the Friedmann equations of cosmology. It 
is easy to show that such cosmological equations follow from Einstein gravity so deriving field equations from a field 
action and then reducing them to cosmological equations or reducing the field Lagrangian to a point -like Lagrangian 
and then deducing the Euler-Lagrange equations gives exactly the same results [12|. In Ea. (l77| . dot represents 
derivative with respect to cosmic time and, as standard for cosmological Lagrangian deduced from field theories, the 
covariance is lost. The variables R and OR can be considered independent and, by the method of Lagrange multipliers, 
we can eliminate higher than one time derivatives. If we would not consider Lagrange multipliers, the Lagrangian 
([77| cannot be considered canonicalflSj. The action related to Lagrangian (I74p . up to OR terms becomes 



A^2tt^ dtl a^F - Ai 



i? + 6 



k 



-A2 



OR-R-'i 



R 



(78) 



A1.2 are given by varying the action with respect to R and OR. that is 
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Only in this case the action results canonically defined in terms of R and OR considered as independent variables. 
After an integration by parts, the point-like Lagrangian results 



. r. ■2dF ^ d fdF\ ,^d f dF \ 
L = baa 77— + 6a a— 77— — a R— „ ,^ 

dR dt \dRj dt \d{DR) J 



^-[^+^)dR-°^dm) 



(80) 



where, clearly, the canonically conjugate variable of configuration space are the set Q — {a, R, DR} and the relative 
velocities. A remark is necessary at this point. We can also take into account 



Ai = a^ 



dF dF 
OR ^ diORj 



(81) 



as a Lagrange multiplier . The Lagrangian which comes out differs from (|80p just for a term vanishing on the 
constraint, being 



£^£~a^R + 6 



□ 



dF 
d{aR) 



(82) 



From this point of view, considering the point-like Lagrangian £ or £ is completely equivalent. 

It is important to stress that Lagrange multipliers are constraints that, after variation, give rise to further equations 
of motion (one for any multipliers). In fact, the expression that are multiplied by Lagrange multipliers in the action are 
constraints. The variation of the action with respect to the Lagrange multipliers gives equations of motion of the form 
"constraint equal to zero". Alternatively, one can solve the Lagrange multipliers and insert them into the action. From 
the resulting action one obtains, of course, the same constraint equations. To show this point, let us derive the Euler- 
Lagrangian equations from the Lagrangian ([50]) that is defined on the tangent bundle TQ = {a, d, R, R, DR, (DR)}. 
They can be also deduced from the Einstein equations (f75|) . The equation for the variables {a, a} gives 



R 



dF 



DR 



2R 



dR 
DR- HR 
d^F 



dF 
d{DR) 
cfF_ 
dR^ ^ 
2 



F 



2{aR) 



R{DR) 
d^F 



3H^ 

d^F 
diORf 



2H 



dF 
dR 



2a^R-2H{DR) 



if 



d^F 



dRdiDR) 



dRd{DRf 



m{DR) 



d^F 



dR^d{DR) 



= 0. 



The equation for {R, (R)} gives 



□ - 



dF 



0. 



d{DR) 

Finally, the equation for the pair {DR, DR} coincides with the Lagrange multipliers 

k 



DR = R + 3HR , 



R 



Q\H + 2H^ 



The energy condition, that is the (0, 0)-Einstein equation, gives 



dF 
dR 



dfdF 



dt \dR 



6 



F-ki 



dF 



dt \d{DR) 



(83) 



(84) 



(85) 



(86) 



This derivation cleary show that dynamics can be made canonical by Lagrange multipliers. However, considering 
further higher-order D^'R terms the process can be made iterative since for each 2-orders (i.e. □) one has another 
Lagrange multiplier. This method allows to select suitable changes of variables that can be identified once the variables 
a, i? and DR are disentangled [12]. As a consequence, dynamics can be reduced and exactly integrated [14, 15]. 

An important remark is in order at this point. If we take into account a conformal transformation as pil) . it is 
easy to show that 



dF\ 
^''"^ I dR) 



(87) 
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for a F(i?)-gravity and 




dUR 



dF 



) 




dUR 



dF 



) 



(88) 



for F{R, □i?)-gravity. It is easy to see that such transformations are related to the Lagrange multiphers ((79|) and 
(I8ip . This means that operating a conformal transformation from the Jordan frame to the Einstein frame, that is 
disentangling the additional gravitational degrees of freedom related to higher order theories of gravity has the same 
dynamical meaning of reducing the dynamical system by imposing Lagrange multipliers. 



In summary, we studied the role of Lagrange multiplier constraint for DE cosmology in scalar-tensor and quintessence 
theory and modified _F'(_R)-gravity. It is demonstrated that the presence of constraint significally simplifies the re- 
construction scenario. Moreover, the details of cosmological evolution are qualitatively changed. For instance, the 
phantom/non-phantom transition in scalar-tensor theory is more smooth. In modified gravity, the presence of La- 
grange multiplier induces the necessity to introduce second F{R) function which plays major role for cosmology. 
Again, the reconstruction program is qualitatively simplified. It is shown that viable reconstruction may be achieved. 
For instance, the examples of dark energy era and unified early-time inflation with late-time acceleration are worked 
out. We considered the constraint of specific form which permits to change the EoS parameter of the effective scalar 
fluid. However, many other choices for constraint may be useful. For instance, additional modification of Lagrange 
multiplier constraint may help to pass the local tests for the theory which cannot pass local tests in its original 
formulation. 

A general comment is in order at this point. Gravitational theories are constrained theories. Such constraints can 
be anholonomic and then can result as further equations of motion for the related dynamical systems. We have flrst 
investigated the possibility that introducing by hand Lagrange multipliers, we can select suitable forms of the effective 
potential V{(j)) and of F(i?)-gravity. This approach leads to solvable DE dynamics for several physically interesting 
models. Besides, we have considered higher-order models. Lagrange multipliers allow, in this case, to make the theory 
canonical, that is they allow to disentangle the degrees of freedom of the problem and then reduce the dynamics. In 
some examples that we have worked out, it is possible to achieve exact solutions thanks to the multipliers that i) 
result as new equations of motion ii) allow suitable change of variables which, identifying cyclic variables, are related 
to constant of motion and allow to reduce and integrate the dynamics. The two methods (i.e. imposing the multipliers 
a priori or using them to reduce dynamics and make it canonical) are effectively equivalent and show the possibility 
to develop a new approach to alternative theories of gravity. 

An important consideration is in order for conformal transformations. In the Einstein frame, gravitational degrees 
of freedom and scalar fleld are well separated. Imposing Lagrange multipliers means to search for suitable forms of 
the scalar field potential that allow to integrate dynamics. In this case, the Lagrange multipliers act as a "selection 
rule" on the dynamics and give fixed stationary points. On the other hand, in the Jordan frame, dynamics is not 
"canonical" since gravitational degrees of freedom and/or scalar fields are not disentangled. Lagrange multipliers, as 
said, make dynamics canonical. It is interesting to see that the form of conformal transformation can be related to 
the Lagrange multipliers (they have the same forms, see Eqs. ([75]) . ([5T|) . ([57)l . and ([55)1 V In some sense we can state 
that imposing Lagrange multipliers in the Jordan frame and performing conformal transformations to the Einstein 
frame are the same operation. Also the reverse statement holds: given a Lagrangian in the Einstein frame endowed 
with Lagrange multipliers means, under conformal transformation to the Jordan frame, to take into account higher 
order or non-minimally coupled theories that are, in any case, canonical. This question will be investigated elsewhere. 
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